Spin-1/2 Optical Lattice Clock 
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We experimentally investigate an optical clock based on ^'^^ Yb (/ = 1/2) atoms confined in an opti- 
cal lattice. We have evaluated all known frequency shifts to the clock transition, including a density- 
dependent collision shift, with a fractional uncertainty of 3.4 x 10~^^, limited principally by uncer- 
tainty in the blackbody radiation Stark shift. We measured the absolute clock transition frequency 
relative to the NIST-Fl Cs fountain clock and find the frequency to be 518 295 836 590 865.2(0.7) Hz. 

PACS numbers: 32.30.-r, 06.30.Ft, 32.70.Jz, 37.10. Jk 
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Optical frequency standards based on the ^5*0 -^^Po 
transition in alkaline earth- (like) atoms confined in an op- 
tical lattice [H have recently demonstrated performances 
approaching or exceeding those of the best microwave 
primary frequency standards Although these 

clocks have not yet achieved the same level of uncertainty 
as the best single trapped-ion clocks they have ad- 
ditional promise in that the large number of atoms si- 
multaneously interrogated (typically tens of thousands) 
should produce lower quantum projection noise, even- 
tually leading to increased clock stability 0, Early 
work on lattice clocks has principally focused on stron- 
tium [H, 0, 01) 9,nd its only abundant fermionic isotope, 
^^Sr, has nuclear spin / — 9/2. Here we present an exper- 
imental study of a lattice clock employing the simplest 
fermionic structure available, a spin-1/2 system. Such a 
system also lends itself well to other investi gati ons, in 
eluding quantum information applications 
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Our previous work on Yb lattice clocks employed a 
bosonic isotope, ^^^Yb (/ — 0) 0, [ll|, which offers the 
simplicity of an effective two-level system. A spin-1/2 
system has only slightly more complicated level struc- 
ture, and reduces the need for careful calibration of 
field-induced shifts that may limit the accuracy of even 
isotope-based clocks Furthermore, a spin-1/2 sys- 

tem offers two advantages over higher angular momen- 
tum systems: optical pumping of atoms to the desired 
771 p ground state is more straightforward, and the tensor 
component of the lattice light shift is absent Al- 
though any clock based on a fermionic isotope needs to 
have its first-order Zeeman sensitivity removed, this is 
readily accomplished by averaging over two transitions 
with opposite dependence (see Fig. la) (1, 13,[3|- 
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Benefitting from the simplicity of the spin-1/2 system, 
our first realization of a ^^^Yb lattice clock has demon- 
strated results competitive with those of the best lat- 
tice clocks to date. By comparing the Yb clock with 
other optical and microwave frequency standards via a 
femtosecond- laser comb 1^ 3] , we have canceled or cal- 
ibrated all known systematic effects with a total uncer- 
tainty of 3.4 X 10~^^ and measured the absolute fre- 
quency with an uncertainty of 1.4 x 10^^^. We have 
observed a density-dependent collision shift that persists 
despite the use of ultracold, spin-polarized fermions. Ad- 
ditionally, in anticipation of higher-dimensional lattice 
geometries that will suppress non-zero density shifts, we 
have measured the vector polarizability of the clock tran- 
sition to be ay « —0.08 a.u. 

The experimental sequence is similar to that we have 
described previously [3^ [HI. ^^^Yb atoms are trapped 
and cooled, first on the strong ^5*0 ^^Pi transition at 
399 nm, and second on the narrower ^5*0 ^ ^Pi tran- 
sition at 556 nm. The atoms are then loaded into a 
one-dimensional, horizontally oriented optical lattice of 
500 Er depth (recoil energy Er/ks = 100 nK). Approxi- 
mately 3 X 10^ atoms with a temperature of 15 /iK are 
captured in several hundred lattice sites, leading to an es- 
timated density of po = 10^^/cm'^. To reduce line-pulling 
effects from other transitions, and in hopes of reduc- 
in g de nsity shifts through the Pauli exclusion principle 
, we pump the population to a single ground state, 
alternating experimental cycles between mp = il/2, be- 
fore exciting the clock transition. This state polarization 
is accomplished by optically pumping on the ^Sq <-^'^Pi 
transition at 556 nm, with which we transfer the atoms 
to a single ground state with 99 % efficiency (Fig. lb). 

During spectroscopy of the ^Sq ^Pq clock tran- 
sition at 578 nm, we apply a small magnetic field of 
0.11 mT. The lattice polarization and probe polariza- 
tion are aligned with the vertical magnetic field, and a 
TT-pulse resonant with the clock transition is applied for 
80 ms, yielding a Fourier-limited lineshape of 10 Hz. The 
578 nm probe light is pre-stabilized to a high finesse, en- 
vironmentally isolated cavity, which produces a fractional 
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FIG. 1: (a) Measured —*^Po spectrum with an externally 
applied magnetic field of B=0.11 mT and 7r-polarized probe 
light. The first-order Zeeman effect splits the two resonances 
by 2f0 Hz/G from center, (b) Zoomed in scan of one of the vr 
transitions, with atoms pumped to (squares) and away from 
(circles) the mp = —1/2 spin state. A Lorentzian fit (solid 
line) gives a Fourier-limited FWHM of 5.5 Ifz (for 160 ms 
probe time), which represents a quality factor of 10^'', with 
over 60 % contrast. Both (a) and (b) were averaged five times 
to suppress clock laser frequency noise. 



instability of « 2 x 10~^^ in 1 s with the drift held be- 
low 1 Hz/s. We detect the fraction of excited atoms by 
applying a sequence of three pulses (5 ms each) resonant 
with the blue cooling transition at 399 nm and collect- 
ing the scattered fluorescence on a photomultiplier tube 
(PMT). Each PMT signal is integrated and then digi- 
tized by a microprocessor. The first pulse measures and 
heats away all atoms remaining in the ground state after 
spectroscopy; the second pulse measures the background 
signal due to fluorescence from the thermal atomic beam; 
and the third pulse measures the atoms that were excited 
to the long-lived '^Pq state by pumping the atoms back 
to the ground state with the ^Pq ^Di transition at 
1388 nm. From these three signals (ground state popula- 
tion, excited state population, background level) the mi- 
croprocessor computes a normalized excitation fraction 
to reduce the effects of shot-to-shot atom number fluctu- 
ations. We alternately probe the half-maximum points of 
the resonance from each mp ~ ±1/2 ground state, and 
the microprocessor controls two synthesizers that lock to 
the average and splitting of the two resonances. This 
detection system is nearly projection noise-limited, with 
a signal-to-noise ratio exceeding 100; however, the clock 
instability is limited by the clock laser 

To evaluate systematic shifts of the clock frequency, we 
performed a series of optical comparisons against the Ca 
optical standard at NIST 18| via a self-referenced fem- 
tosecond laser frequency comb [isl . [lij . During these 
measurements, the fractional frequency instability be- 



tween the two clocks started at 5 x 10""'^^ at 1 s and 
averaged down to 3 x 10^^^ at 200 s. On longer time 
scales, drifts in the Ca clock degrade the stability, so typ- 
ically data sets were taken in repeated bins of 200 s. 

As with all optical lattice clocks, it is necessary to mea- 
sure the effects of the Stark shift from the confining lat- 
tice light. The flrst-order effect is canceled by tuning the 
lattice to the "magic" frequency, at which the ground and 
excited clock states have identical polarizability values 
[l[. To flnd this magic frequency, we varied the lattice 
intensity to find the Stark shift for a given lattice fre- 
quency, and then repeated this procedure for a series of 
frequencies. The results of this measurement are shown 
in Fig. [2^, and the magic frequency is determined to be 
394 798 329(10) MHz. This gives a difference between 
magic frequencies for the two measured isotopes to be 
^}^a^^c - ^'rIag^c = 1146(36) MHz 19]. Wc also uotc that 
this measurement is complicated by the density shift in- 
duced by a change in the lattice confinement. The final 
uncertainty for the scalar Stark shift is given by the mea- 
surement statistics at the magic frequency (1.7 x 10""'^^), 
which could be substantially reduced by measuring with 
higher lattice power (20| to achieve a larger lever arm, 
with an additional term (1 x 10~^^) allowing for possi- 



ble magnetic dipole and electric quadrupole shifts [21 1 



In addition to the first-order scalar Stark shift, there 
are two other lattice effects that need to be considered. 
First, two-photon resonances near the magic frequency 
lead to hyperpolarizability shifts; these were addressed 
in [l9|. The shift and uncertainty for our typical lattice 
depth are shown in Table HI Second, for odd isotopes, 
vector and tensor light shifts arise from hyperfine inter- 
actions, though for spin-1/2 systems the tensor shift is 
zero [13]. As expressed in [22|], the vector shift for a given 

mp state is Sfy — —mpAay (^fo) , where A is the de- 
gree of circular polarization, Eq the electric field, and a„ 
the vector polarizability at 759 nm. For our experimental 
geometry, the vector and Zeeman shifts of the ^Pq add in 
quadrature, allowing the vector shift to be determined by 
measuring the splitting between mp = ±1/2 resonances 
as a function of lattice polarization ellipticity. As shown 
in Fig. the maximum vector shift (for circular po- 
larization with 50 /iK lattice depth) is Sf^ =220(30) Hz, 
where the error bar is derived from imperfect polarization 
control due to birefringence in the vacuum windows. This 
shift corresponds to a vector polarizability of w —0.08 
a.u., in agreement with the calculation in (2^ . During 
typical operation, the vector shift is suppressed to less 
than 10~^^ by three effects: linear polarization [A « 0), 
geometry [22|, and cancellation by averaging over both 
mp — ±1/2 resonances. For a three-dimensional lat- 
tice, the vector shift could be more problematic, as each 
lattice site could experience a different, time-dependent 
composite polarization vector. This can lead to broad- 
ening of each resonance by as much as Sfv, requiring the 
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FIG. 2: Measured systematic eflFects. (a) Lattice Stark shift. 
By varying the lattice power during spectroscopy for differ- 
ent lattice wavelengths, we find the magic frequency to be 
394,798,329(10) MHz. A linear fit (solid line) gives the slope 
as -11(1) Hz/GHz. (b) Vector Stark shift. By monitoring the 
splitting between the two tt resonances as the lattice elliptic- 
ity is varied, we find an upper limit for the vector Stark shift 
to be 220(30) Hz for a 50 {xK lattice depth, (c) Density shift 
as a function of final excitation, which was varied by changing 
the lock point on the resonance, (d) Magnetic field. Measur- 
ing the clock frequency versus external magnetic field, we find 
the first-order Zeeman to be canceled at the 4 x 10~^^ level, 
while the second-order coefficient of -0.07(1) Hz/G^ is in good 
agreement with previous work [3, [l^l ■ A quadratic fit (solid 
line) is also shown. 



vector shift to be suppressed by other means (e.g. [23j). 
Even with vector shift-induced broadening, the clock fre- 
quency remains unperturbed by proper averaging over 
both mp = ±1/2 resonances. 

To measure a possible density-dependent collision 
shift, we lowered the power in the 399 nm slowing beam 
to reduce the number of atoms loaded into the lattice 
while looking for a clock shift. Sideband spectra show 
that the atom temperature (« 15 /iK) is independent of 
slowing power, ensuring that only the atom number (not 
the confinement) is varied. We find a linear relation- 
ship between atom number and clock frequency, with a 
slope of -0.85(4) Hz/po- We emphasize that while abso- 
lute density is difficult to measure, the relative density is 
well-calibrated by the detection scheme described above. 
During typical operation, the excitation fraction is 0.3, 
at which the density shift uncertainty is 8xl0~^^. 

We also investigated the density shift dependence on 
atomic excitation by varying the detuning of the lock 
points on the clock transition (see Fig. 2c). Because the 
magnitude of the shift increased as the atomic excita- 



tion was reduced, we first investigated two-body colli- 
sions between ^5*0 atoms as the dominant shift mecha- 
nism. Founded in the Pauli exclusion principle, as two 
ground {^Sq) state atoms in the same spin state comprise 
a pair of identical fermions, s-wave interactions are pro- 
hibited. The atomic temperature is not sufficiently cold 
to rule out by itself the p-wave term of such an inter- 
action. However, photoassociation spectroscopy of many 
Yb isotopes has been used to determine cold collision 
scattering cross sections of ground state atoms, not only 
yielding precise s-wave scattering lengths but also en- 
abling calculation of higher wave contributions 2^ . The 
expected interaction for -'^'^^Yb is not sufficiently large to 
justify the observed shift. Furthermore, the interaction 
is expected to be much larger for ^^''Yb, where no coUi- 
sional shift has been observed at similar densities (albeit 
with larger uncertainty) Q. 

Another possibility lies with interactions between 
and ^Pq atoms. In ultracold fermionic Sr, inhomogeneous 
excitation of the lattice-confined atoms leads to only par- 
tial suppression of s-wave interactions of this type and 
has been shown to yield a non-zero density shift [l3|- 
This same framework, applied to ^^^Yb, including only 
partial suppression by the two body correlation function, 
appears consistent with our observed shift as a function 
of excitation. In contrast to the case of Sr, here no clear 
zero crossing in the shift is present. However, based on 
the above analysis, this zero crossing is expected to oc- 
cur at higher excitation fraction, and the crossing would 
be found at reduced excitation for a colder atomic sam- 
ple. Colder temperatures would also give overall smaller 
shifts. These shifts would be further reduced by sparse 
population in a higher dimensional lattice. 

A summary of all known systematic effects is given 
in Table HI We determined the Zeeman shifts by mea- 
suring changes in the clock frequency as the magnetic 
field is varied. The resulting data, shown in Fig. 2d, 
are fit to linear and quadratic functions to determine the 
first- and second-order sensitivities. As expected, the 
first-order effect is canceled by the averaging technique, 
while the second-order shift is in good agreement with 
theory [l3| and previous measurements 0]. The probe 
light shift caused by off-resonant excitations originating 
from the 578 nm light was calculated based on previous 
measurements [s] for our typical 20 nW of power. The 
acousto-optic modulator (AOM) phase chirp was mea- 
sured by changing the RF power to the chopping AOM 
(between cavity and atoms) and lookin g fo r a clock shift, 
as well as by heterodyne interferometry 25]; both of these 
measurements were consistent with no shift for our typi- 
cal conditions within the measurement uncertainty. The 
blackbody shift for room-temperature operation is calcu- 
lated based upon the polarizability calculations in [2^. 
The 10 % uncertainty in the polarizability leads to a 
room-temperature uncertainty of 2.5 x 10~^^, and future 
work to measure the polarizability or operate below room 
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TABLE I: Uncertainty budget for systematic effects in the 
^'^^Yb optical lattice clock. Values are given in fractional units 
10~^® and represent typical operating conditions. 



Effect 


Shift Uncertainty 


Blackbody 


-25 


2.5 


Lattice polarizability 


0.4 


2.0 


Hyperpolarizability 


3.3 


0.7 


Density 


-16.1 


0.8 


First-order Zeeman 


0.4 


0.4 


Second-order Zeeman 


-1.67 


0.1 


Probe light 


0.05 


0.2 


AOM phase chirp 


-0.05 


0.1 


Servo and Line pulling 





0.1 


Residual Doppler 





0.1 


Total 


-38.7 


3.4 



temperature will be necessary to improve the clock's ac- 
curacy. Alternatively, one could define the fundamental 
clock operating condition at room-temperature, in which 
case the correction is much smaller and the blackbody 
uncertainty is 3.3 x 10~^^ for 1 K uncertainty in tem- 
perature. Other systematic effects (servo error, resid- 
ual Doppler shifts) have been calculated to be negligibly 
small. 

We have completed a comparison of the ^^^Yb lattice 
clock with the NIST-Fl Cs primary standard 0. For 
these measurements, the femtosecond comb's repetition 
rate 1 GHz) was locked to the heterodyne beat with 
Yb, mixed down to several MHz with a synthesizer ref- 
erenced to a hydrogen maser, and counted, while the Cs 
clock simultaneously calibrated the H-maser. Allan devi- 
ation and histogram plots of the measurement statistics 
give the fractional instability as 5.1 x 10~^^t^^^^ . Af- 
ter averaging for > 10^ seconds, we find the measured 
frequency is 518 295 836 590 865.2(0.7) Hz, with the sta- 
tistical uncertainty (1.4 x 10~^^) still much larger than 
the systematic uncertainties of either clock. Combined 
with our previous measurements we determine the 
isotope shiituZik - ^Tck = 1 811 281 647.4(1.1) Hz. 

In conclusion, we have completed a systematic evalu- 
ation of the ^^^Yb optical lattice clock with uncertainty 
3.4 X 10^^®, competitive with the best neutral-atom mi- 
crowave and optical frequency standards currently oper- 
ating. The simplicity of the spin-1/2 system is enjoyed in 
several ways, including the simple spectrum, straightfor- 
ward state polarization, and lack of tensor lattice shifts. 
These results indicate that ^^^Yb is an excellent candi- 
date not only for an optical frequency standard but also 
for quantum information experiments that could bene- 
fit from a spin-1/2 quantum system with long coherence 



times and simple nuclear spin manipulation. 

The authors thank L.-S. Ma for experimental assis- 
tance, S. Bickman and T. Johnson for careful reading of 
the manuscript, and P. Julienne, S. Blatt, and J. Sher- 
man for helpful discussions. 

Note added. — While completing this manuscript, 
we became aware of another recent measurement of the 
^''^Yb clock frequency [2§|, which agrees with our value 
to within their stated uncertainty. 
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